Journal  of  Power  Sources  188  (2009)  475-482 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


pri 

SllbudtS 


Control-oriented  thermal  management  of  solid  oxide  fuel  cells  based  on  a 
modified  Takagi-Sugeno  fuzzy  model 

Jie  Yanga  b,  Xi  Lic  *,  Hong-Gang  Moub,  Li  Jiana 

a  School  of  Materials  Science  and  Engineering,  State  Key  Laboratory  of  Material  Processing  and  Die  &  Mould  Technology, 

Huazhong  University  of  Science  &  Technology,  Wuhan  430074,  China 

b  School  of  Mechanical  and  Electronic  Information,  China  University  of  Geosciences,  Wuhan  430074,  China 

c  Department  of  Control  Science  and  Engineering,  Key  Laboratory  of  Education  Ministry  for  Image  Processing  and  Intelligent  Control,  Huazhong  University  of  Science  &  Technology, 
Wuhan  430074,  China 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  2  December  2008 
Accepted  4  December  2008 
Available  online  9  December  2008 


Keywords: 

Solid  oxide  fuel  cell  (SOFC) 
Takagi-Sugeno  (T-S)  fuzzy  model 
Control-oriented  modeling 
Enthalpy-balance  equation 


Thermal  management  for  a  solid  oxide  fuel  cell  (SOFC)  is  actually  temperature  control,  due  to  the  impor¬ 
tance  of  cell  temperature  for  the  performance  of  an  SOFC.  An  SOFC  stack  is  a  nonlinear  and  multi-variable 
system  which  is  difficult  to  model  by  traditional  methods.  A  modified  Takagi-Sugeno  (T-S)  fuzzy  model 
that  is  suitable  for  nonlinear  systems  is  built  to  model  the  SOFC  stack.  The  model  parameters  are  initialized 
by  the  fuzzy  c-means  clustering  method,  and  learned  using  an  off-line  back-propagation  algorithm.  In 
order  to  obtain  the  training  data  to  identify  the  modified  T-S  model,  a  SOFC  physical  model  via  MAT- 
LAB  is  established.  The  temperature  model  is  the  center  of  the  physical  model  and  is  developed  by 
enthalpy-balance  equations.  It  is  shown  that  the  modified  T-S  fuzzy  model  is  sufficiently  accurate  to 
follow  the  temperature  response  of  the  stack,  and  can  be  conveniently  utilized  to  design  temperature 
control  strategies. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC),  operated  at  temperatures 
between  600  and  1000  °C,  is  a  promising  technology  for  sustainable 
energy  conversion,  having  a  wide  range  of  possible  applications.  Its 
unique  advantages  are  high  efficiency,  fuel  flexibility,  and  low  pol¬ 
lutant  emissions,  while  combined  heat  and  power  generation  can 
be  realized  with  the  waste  heat  produced  in  SOFCs  [1].  However, 
the  high  operating  temperature  also  brings  serious  challenges  to 
the  technology,  such  as  material  selection,  thermal  management, 
reliability  and  durability.  Thermal  management  associated  with 
temperature  and  its  distribution  is  a  critical  factor  for  the  cell  per¬ 
formance  and  stack  behavior,  affecting  cell  integrity,  cell  voltage, 
current  density  distribution  and  power  output  [2,3]. 

An  SOFC  is  a  nonlinear  multi-variable  system  which  is  difficult  to 
model  by  traditional  methodologies.  In  order  to  effectively  control 
the  cell  temperature  in  terms  of  control  strategies,  it  is  necessary 
to  develop  model-based  control  methods.  In  recent  decades,  math¬ 
ematical  models  have  been  established,  emphasizing  the  accurate 
description  of  cell  internal  mechanisms,  such  as  mass-  and  energy- 
balances  and  electrochemical  kinetics  [3-8].  These  models  are 
useful  for  cell  design,  analysis  and  optimization.  However,  they 
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are  complicated  for  use  as  model-based  control  methods  [9],  and 
in  some  cases  they  are  impracticable.  In  addition,  researchers 
have  tried  to  establish  models  via  data-driven  approaches  so  as 
to  obtain  valid  control  strategies  [10,11].  To  control  the  complex 
SOFC  system,  simpler  data-driven  models  based  on  system  identifi¬ 
cation  methodology  are  useful  and  convenient.  Arriagada  presented 
an  artificial  neural  network  to  identify  the  SOFC  model,  which 
was  trained  by  experimental  data  for  prediction  of  various  SOFC 
operational  parameters.  This  method  can  accomplish  performance 
analysis  and  prediction,  but  it  is  not  perfect  for  the  control  purpose, 
and  the  experimental  data  may  not  sufficiently  reflect  all  the  char¬ 
acteristics  of  an  SOFC  [12].  Jurado  [13,14]  proposed  a  Hammerstein 
model  for  identification  and  predictive  control.  This  model  is  specif¬ 
ically  suitable  for  nonlinear  systems  where  the  nonlinear  block  is 
static  and  followed  by  a  linear  system,  and  can  only  be  applied  for 
small  signal  and  transient  stability  studies.  Huo  et  al.  [15]  and  Kang 
et  al.  [16]  used  the  least  squares  support  vector  machine  to  establish 
a  data-driven  model  for  an  SOFC,  achieved  a  solution  by  solving  a  set 
of  linear  equations.  This  model,  based  on  mass-  and  energy-balance 
equations,  describes  the  temperature  dynamics  and  predicts  stack 
voltage. 

Fuzzy  modeling  and  identification  based  on  measured  data 
are  effective  tools  for  the  approximation  of  an  uncertain  non¬ 
linear  system  [17].  Among  various  fuzzy  modeling  methods,  the 
Takagi-Sugeno  (T-S)  fuzzy  model  has  attracted  most  attention 
[18]  due  to  it  having  fewer  rules  than  other  models.  Each  rule 
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Nomenclature 

A  cell  active  area  (m2) 

Cp  specific  heat  capacity  of  gas  species  (J  mol-1 1<-1 ) 
Cp  area  specific  heat  capacity  (J  K-1  m-2) 

F  Faraday’s  constant  (96,485  C  mol-1 ) 

h  specific  enthalpy  (J  mol-1 ) 

H  enthalpy  (J) 

H  enthalpy  flow  (W) 

i  current  density  (A  m-2 ) 

i0  exchange  current  density  (A  m-2 ) 

ir  reaction  current  density  (A  m-2 ) 

I  current  (A) 

m  mass  (kg) 

N  molar  flux  (mol  m-2  s-1) 

n  number  of  moles  (mol) 

h  molar  flow  (mol  s-1) 

rii  molar  number  of  species  i  (mol) 

hc  hydrogen  combustion  molar  flow  (mol  s-1 ) 

P  pressure  (bar) 

R  universal  gas  constant  (8.314 J  mol-2  K-2 ) 

t  time  (s) 

T  temperature  (I<)  also  period  (s) 

T0  ambient  temperature  (I<) 

U  voltage  (V) 

V  volume  (m2) 

x  input  vector 

Xi  molar  fraction  of  species 

y  output  vector 

Greek  letters 

a  symmetry  factor  learning  factor 

a  T-S  factor 

p  overpotential  (V) 

v  stoichiometric  coefficient 

Subscripts 

i  gas  species 

r  reaction 

Superscripts 
in  fuel  cell  inlet 

out  fuel  cell  outlet 


consequence  with  a  linear  function  can  be  described  by 
input-output  mapping  in  a  wide  range  and  the  fuzzy  implication 
used  in  the  model  is  relatively  simple  [19].  The  T-S  fuzzy  model  has 
the  capability  to  approximate  a  large  class  of  static  and  dynamic 
nonlinear  systems  [20,21],  and  has  been  successfully  proved 
to  be  applicable  in  model  identifications  and  control  strategies 
[22-24].  In  most  previous  applications,  single-input/single-output 
or  multi-input/single-output  was  considered  [25-27].  Yang  et  al. 
[28]  described  a  dynamic  T-S  fuzzy  model  of  a  molten  carbonate 
fuel  cell,  consisting  of  voltage  and  temperature  submodels,  but  the 
temperature  submodel  can  only  be  valid  within  a  relatively  narrow 
temperature  range.  The  modeling  of  SOFCs  using  the  T-S  model  has 
not  been  reported  so  far. 

In  this  study,  a  control-oriented  thermal  management  of  an 
SOFC  stack  using  a  modified  T-S  fuzzy  model  was  conducted  with 
an  identification  method.  The  training  data  were  extracted  from 
an  SOFC  stack  physical  model  that  was  constructed  with  MATLAB. 
The  physical  model  built  by  mass-  and  energy-balance  equations, 
and  enthalpy-balance  equations  can  serve  as  a  real  SOFC  stack.  The 
temperature  subsystem  based  on  enthalpy-balance  equations  can 


realize  temperature  control  strategy.  This  T-S  fuzzy  model  consists 
of  if-then  rules  with  fuzzy  antecedents  and  mathematical  func¬ 
tions  in  the  consequent  part.  Unlike  the  original  T-S  model,  SOFC 
input  and  output  variable  spaces  were  divided  into  many  subspaces 
according  to  the  principle  of  the  fuzzy  c-means  (FCM)  clustering 
method  [29].  The  model  parameters  were  learned  by  using  the 
back-propagation  (BP)  learning  algorithm  [30].  The  task  of  system 
identification  was  to  determine  both  the  nonlinear  parameters  of 
antecedents  and  the  linear  parameters  of  the  rule  consequence.  The 
satisfactory  results  show  that  the  modified  T-S  fuzzy  model  can 
achieve  the  temperature  control  strategy. 

2.  SOFC  physical  model 

In  this  SOFC  physical  model,  some  simplifications  and  assump¬ 
tions  are  made,  since  high  accuracy  is  not  necessary  for  a  physical 
model  that  serves  for  a  subsequent  control  strategy.  Any  devia¬ 
tions  between  the  model  and  the  real  fuel  cell  can  be  managed 
by  a  feedback  loop  in  the  control  system  [31].  The  proposed  SOFC 
stack  physical  model  has  the  following  assumptions: 

(1 )  Stack  is  fed  with  hydrogen  and  air;  the  fuel  processor  dynamics 
are  not  included. 

(2)  A  uniform  gas  distribution  among  the  cells  is  assumed,  since 
there  is  a  small  deviation  of  the  gas  distribution. 

(3)  There  is  no  heat  transfer  among  the  cells.  Each  cell  has  the  same 
temperature  and  current  density. 

(4)  There  is  no  heat  exchange  between  the  stack  and  the  ambient 
environment. 

(5)  The  channels  that  transport  gases  along  the  electrodes  have  a 
fixed  volume  and  small  dimension,  so  there  is  a  constant  pres¬ 
sure  in  the  stack. 

2.1.  Mass-balance  model 


For  a  generic  species  i,  the  dynamic  mole-balance  is: 


_  zAnJin  {.out 

dF-n  x>  _n 


(1) 


The  stoichiometric  factor  indicates  how  many  moles  of  the 
species  are  produced  or  consumed  for  each  mole  of  electrons 
transferred.  The  anodic  and  cathodic  reactions  in  the  SOFC  are, 
respectively: 


H2  +  02“^  H20  +  2e“ 


(l/2)02  +  2e-^  O2-  (2) 

then  vq2  =  -V2-  vh2  =  -1,  vHl0  =  1,  vn2  =  0. 

According  to  the  Butler-Volmer  equation  [32]: 

jr  =  fo(ea(nF/RT)r]  _  g-(l-a) (nF/RT)rjj  (3 ) 


When  pressure  is  constant,  the  total  outlet  molar  flow  depends 
only  on  the  transient  reaction  rate  and  temperature  and  can  be 
calculated  using  the  following  relation: 


■out  -in  FA  pV  dT 

nout  =  n,„  +  ^_  +  __ 


(4) 


The  term  (irA/F)  represents  the  algebraic  sum  of  gas  reac¬ 


tion  molar  flow,  and  ( pV/RT2)(dTldt )  reflects  thermal  expansion. 
Inserting  Eq.  (4)  into  Eq.  (1)  and  assuming  constant  pressure,  the 
molar-balance  becomes: 


pV  d.X[ 
RJ~dt 


(x'n  -  Xj)  +  ^(Vj  -  *£>) 


(5) 
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2.2.  Energy -balance  model 

Energy-balance  is  a  very  important  aspect  for  SOFCs  operated 
at  high  temperature.  Here  only  the  enthalpy-balance  is  considered. 
The  main  sources  and  sinks  of  heat  for  an  SOFC  are  entering  and 
exiting  flow,  the  heat  generated  by  cell  reactions,  and  the  heat  lost 
to  the  environment.  An  enthalpy-balance  equation  yields: 


d  lj 

=  H"1  -  Hout  -  irVA  -  H1 
at 


(rloSS 


(6) 


The  term  irVA  represents  the  electrical  power  generated.  The  rela¬ 
tionship  between  enthalpy  and  temperature  can  be  expressed  as 


dH  A.  dT 
dt  p  dt 


(7) 


The  entering  enthalpy  flow  also  contains  the  heat  generated  by 
the  combustion  of  hydrogen,  hence  the  total  entering  enthalpy  flow 
is  equal  to  the  enthalpy  associated  with  the  ambient  air  plus  the 
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Fig.  1.  MATLAB  diagram  of  air  and  H2  flowing  into  the  Stack  Subsystem  via  the  supply  subsystem. 
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entering  hydrogen  flow: 


Hin  =  hairhair(T0)  +  (hn2  +  hcH2)hH2(TH2)  (8) 

In  reality,  Th2  may  be  significantly  different  from  T0,  as  hydro¬ 
gen  from  various  sources  may  undergo  different  processes,  such 
as  expansion,  evaporation  or  dissociation;  however,  for  simpli¬ 
fication  Th2  =  To  is  assumed  without  causing  any  meaningful 
deviation. 

The  total  outlet  molar  flow  will  be  larger  than  or  equal  to  the 
entering  one  and  the  expression  of  this  is: 

Hout  =  ^ri°uth,(T)  (9) 


and  the  outlet  molar  flow  for  species  i,  h°ut  can  also  be  expressed 
as 


nuui  =  nm  +  v 


("&2  +  +*i 


pV  dT 
Kp~dt 


(10) 


where  h|.n  is  the  entering  molar  flow,  +  (irA/2F))  is  the  molar 
flow  associated  with  hydrogen  reaction  and  combustion,  and  the 
term  ( pV/RT2)(dTldt )  is  related  to  gas  thermal  expansion. 

The  heat  loss  to  the  environment,  Hlosst  is  generally  a  benefi¬ 
cial  factor,  as  a  fuel  cell  stack  prefers  heat  to  be  removed  from 
it.  However,  its  actual  expression  depends  on  the  stack  materials 
and  layout.  For  the  sake  of  simplicity,  it  is  assumed  that  Hloss  =  0 
[31]. 
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Fig.  2.  MATLAB  diagram  showing  air  and  H2  flow  from  the  ca_manifold  and  an_manifold  into  the  temperature  block. 
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Fig.  3.  The  structure  of  the  temperature  model  subsystem. 


2.3.  Temperature  control  dynamic  model 

The  objective  of  temperature  control  is  to  keep  the  stack  at  a 
given  temperature  by  manipulating  the  air  flow.  Having  assumed 
that  temperature  is  uniform  in  a  stack,  Eq.  (6)  determines  its  dynam¬ 
ics.  According  to  simultaneous  Eqs.  (6)  and  (7),  it  can  be  obtained 
that: 

Acff-  =  Hin  -  Hout  -  irVA 

dt 

=  (Hin(T)  —  H0Ut(T))  -  (Hin(T)  —  Hin(To))  -  irVA  (11) 

The  terms  {Hm{T)  -  H0Ut{T))  and  {Hin{T)  -  Hm{T0)),  respectively, 
represent  the  reaction  heat  and  the  sensible  heat. 

In  Eq.  (11)  it  assumed  that  the  specific  heat  capacity  for  all  the 
species  is  approximately  the  same  on  condition  that  the  pressure 
is  constant.  This  is  adequate  for  the  control-oriented  modeling. 
According  to  Eq.  (11),  the  temperature  control  model  can  be  built 
with  MATLAB. 

2.4.  Temperature  control  model  realization  in  MATLAB 

The  stack  temperature  is  influenced  essentially  by  the  entering 
flows  and  the  reaction  rate,  and  the  reaction  rate  is  determined 
according  to  the  requirements  of  the  electric  utility.  Therefore,  the 
only  remaining  way  to  influence  the  stack’s  temperature  is  chang¬ 
ing  the  entering  flows.  This  is  the  standard  operating  mode,  where 
the  entering  flows  are  manipulated  to  approach  the  stack  reference 
temperature.  This  physical  MATLAB  model  has  two  input  variables, 
air  flow  and  H2  flow,  and  one  output  variable,  T.  Once  more  vari¬ 
ables  are  selected  from  the  physical  model  and  more  training  data 
are  obtained  for  the  modified  T-S  model,  the  multi-input  (load 
current,  air  flow  hydrogen  flow  rate,  etc.)  and  multi-output  (stack 
temperature,  stack  voltage  stack  power,  etc.)  system  can  be  estab¬ 
lished. 

A  virtual  3.5  kW  SOFC  stack  is  used  in  the  simulation.  The  stack 
consists  of  60  cells  with  anode  and  cathode  gases  in  cross-flow 
and  cell  active  area  0.01  m2.  As  shown  in  Figs.  1-3,  the  physical 
model  replaces  the  real  SOFC  stack  to  generate  the  simulation  data 
required  for  the  modified  T-S  fuzzy  model.  In  Fig.  1,  the  input  vari¬ 


ables  are  air  flow  and  H2  flow.  Air  and  H2 ,  respectively,  enter  into  the 
Air  Supply  Subsystem  block  and  the  Fuel  Supply  Subsystem  block. 
The  internal  structure  of  the  Stack  Subsystem  block  is  shown  in 
Fig.  2.  The  Stack  Subsystem  block  includes  a  cathode  flow  block, 
anode  flow  block,  stack  block  and  temperature  block.  Eq.  ( 11 )  is  built 
in  the  temperature  block  shown  in  Fig.  3;  the  stack  temperature 
(Temp .stack)  can  be  obtained  according  to  Eq.  (11). 

The  enthalpy  expressions  can  be  acquired  from  the  published 
literature  [33]  and  the  relationship  between  enthalpy  and  tem¬ 
perature  for  different  gases  is  reflected  in  these  expressions.  The 
Enthalpy  .out  block  and  Enthalpyjn  block  can  be  expressed  as  fol¬ 
lows  (w  expresses  molar  flow): 

•  For  Enthalpyjn  block: 

h An  =  w.air  x  ft.air  +  w_H2  x  h_H2  +  w.V  x  fr.V  (12) 

where 

ft  .air  =  -1 .0947e4  +  32.50  x  T_air 

ft_H2  =  -0.9959e4  +  30.73  x  T.fuel  (13) 

ft.V  =  -25.790e4  +  42.47  x  T_fuel 

TT  w_fuel  x/Lmoelra  x  2.016 

w.H?  =  - - - 

(fi- moelra  x  2.016  +  18.02)  (14) 

w.V  =  w.fuel  -  w_H2 

•  For  the  Enthalpy  .out  block: 

fr.out  =  w_02  x  h_02  +  w_N2  x  h_N2  +  w_H2  x  h_H2 

+w_V  x  h.V  (15) 


and 

h_N2  = 

— 1 .0590e4  +  31 .40  x  T 

h.  H2  = 

-0.9959e4  +  30.73  x  T 

hj 02  = 

-1.2290e4  + 35.12  x  T 

h.V  = 

-25.790e4  + 42.47  x  T 
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Fig.  4.  The  temperature  response  of  the  SOFC  stack  physical  model. 


Based  on  the  physical  model,  it  can  be  established  that  the  load 
signal  change  can  vary  the  SOFC  stack  temperature.  In  addition,  the 
stack  temperature  also  fluctuates  with  the  input  variables  (air  flow 
and  H2  flow).  As  shown  in  Fig.  4,  the  temperature  of  the  SOFC  stack 
fluctuates  between  815  and  925  K  with  an  ascending  trend  and  the 
temperature  curve  is  fluctuating  rather  than  smooth.  Sharp  fluc¬ 
tuations  of  stack  temperature  present  difficulties  for  temperature 
control,  and  therefore  a  data-driven  model,  instead  of  a  physical 
model,  is  used  to  carry  out  the  control  strategy.  For  this  purpose, 
a  modified  T-S  model  is  applied  with  data  provided  by  the  above 
physical  model. 

3.  A  modified  T-S  fuzzy  model  of  the  structure  and 
algorithm 

3.1  The  structure  of  T-S  fuzzy  model 

An  SOFC  system  can  be  considered  as  a  multi-input/multi¬ 
output  system  that  can  be  divided  into  q  (number  of  input  variables) 
multi-input  and  single-output  (MISO)  systems,  and  each  MISO 
system  can  be  fitted  by  a  fuzzy  T-S  model  independently.  The 
discrete-time  system  can  be  modeled  by  the  collection  of  c  fuzzy 
rules.  It  has  been  theoretically  proved  that  the  T-S  model  has  the 
capability  of  closely  approximating  a  real  system  by  increasing  the 
number  of  fuzzy  rules;  therefore,  the  model  is  highly  precise. 

The  modified  T-S  model  presented  here  is  a  fuzzy  logic  system 
which  is  jointly  composed  of  the  fuzzy  rule  base  and  weight  mean 
defuzzier. 

The  fuzzy  rule  is  in  the  form  of  “if. .  .then. . .”.  The  ith  rule  of  the 
Ith  output  yifk  + 1 )  is  given  by 

Ru:Ifx{k)isAlp  then 

ni(^+^)  =  P-0+P-l^l  +---+P-n*kn  (i  =  1,  •  •  •,  c)  V  ' 

where  c  is  the  number  of  rules,  Al  =  {a!  v  is  the  set 

of  membership  functions  associated  with  the  ith  rule,  and  pi  = 
[pi  0,  p\  1 ,  •  •  •,  p\  n\  is  the  parameter  vector  of  the  ith  submodel.  It  is 
defined  that: 

{UJ  =  [ ui  (k)  >  •  •  •>  uj  (k  -  nuj  +  1 )]  (/=!>•• nj) 

{y;(k)l  =  [y;(k),  ■  ■  ■,  yi(k  -  nyl  +  1 )]  (/  =  1 ,  ••  •,  (?) 

nuj  and  nyl  are  the  order  of  the  input  variable  Uj  and  output  variable 
yh  respectively,  and  each  MISO  subsystem  can  be  denoted  as 


(18) 


Wherex(k)=  { [m  (fc)]gul ,  ■  ■ [unj\™j ,  [yi(k)]^1,  [y,(k)]™}  = 

[Xfei,  ■  ••jXfen]  is  the  regression  data  vector  consisting  of 
input-output  data  at  the  kth  instant  and  before. 

The  antecedent  part  of  “If. . .”  in  the  form  of  the  fuzzy  set  is 
equal  to  the  fuzzy  partition  of  data  space  of  the  SOFC  input-output 
variables,  and  the  consequent  part  of  “then. . .”  represented  by  a  lin¬ 
ear  functional  relation  is  a  linear  composition  of  the  input-output 
variables.  The  identification  is  to  find  the  fuzzy  SOFC  model  with 
reasonable  complexity  to  minimize  the  error  between  the  fuzzy 
SOFC  model  output  and  the  real  SOFC  stack  output.  In  this  model, 
what  is  to  be  partitioned  is  the  whole  space  spanned  by  regression 
data  vectors,  but  not  the  single  variable  space  presented  by  Takagi 
and  Sugeno  [18]. 

Membership  function  of  the  fuzzy  set  Aj  can  be  defined  as  the 
product  of  Gaussian  functions: 


where  xm  is  the  mth  regression  data  vector,  p(=l, . .  .,n)  is  the  dimen¬ 
sion  of  data  vector,  V\ =[v^t  vi2 ,  . . .,  vin\  and  crj=[orn,  crl2 ,  . . crin] 
are  subsequent  identification  parameters  of  the  membership  func¬ 
tions.  According  to  the  above  simultaneous  equations,  the  fuzzy 
reasoning  output  can  be  expressed  as 


yi(k  +  i) 


with 


ELiw.y;,i(k  +  1) 

Ei=iw< 

Et  1  W.'(p!  0  +  p\  V *kl  +  ■  • '  +  p'i  nXk "  ) 

EL,wf 

EtlW'(U  xml  xp[)  _  fa 

Eli  wi  t2 


1  f  X-mp  —  Vjp 

2  V  ajp 


2 


(20) 


(21) 


where  [p\  0-  ■  plc  0-  •  p\  r  •  plcX  ■  p\  n-  •  -p1^  is  the  subsequent 
identification  parameter  vectors  of  the  consequent  part,  and  wz  is 
the  match  degree  of  mth  data  vector  corresponding  to  ith  rules. 


3.2.  Model  parameters  of  BP  learning  algorithm 

The  above  fuzzy  logic  system  can  be  represented  as  a  feed¬ 
forward  network  in  Fig.  5,  whose  weight  coefficients  have  clear 
physical  meanings  and  can  be  learned  by  back-propagation  of  the 
output  predictive  error  [23,30].  A  parameter  set  of  the  fuzzy  system 
that  minimizes  the  output  predictive  error  can  be  found  through  a 
collection  of  training  data  produced  by  the  physical  SOFC  model: 

£ ^—ym=l m  yi,m)  (22) 

where  y/  m  is  the  SOFC  physical  model  output  of  mth  data  vector,  and 
yi  m  is  predictive  output  of  fuzzy  inference  system.  Calculating  the 
gradient  of  predictive  output  error  E  through  Eqs.  (21 )  and  (22)  can 
be  used  to  obtain  the  learning  algorithm  of  fuzzy  inference  system 
parameter  vectors  uz,  az  and  pp. 

vUp{k+\)=vUp{k)-aff-{k)  (23) 

3£ 

aUp(k+\)  =  aUp(k)-a-^-(k) 


y,(k+l)=/[x(k)]  l  =  1,  -,q 


(24) 
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1  Xmn 

Fig.  5.  Feed-forward  network  form  of  fuzzy  logic  system. 


Table  1 

Parameters  of  the  SOFC  stack  used  in  the  physical  model. 


Item 

Value 

Cell  number 

60 

Cell  active  area  (m2) 

0.01 

The  heat  capacity  per  cell  area  (J  I<-1  m-2) 

7000 

The  arm  power  (kW) 

3.5 

Anode  total  volume  (m3) 

0.005 

Cathode  total  volume  (m3) 

0.001 

Gas  inlet  flow  constant 

0.3629  x  10“5 

Anode  gas  outlet  flow  constant 

0.0544  x  10“5 

Cathode  gas  outlet  flow  constant 

0.2118  x  10“5 

Air  heat  capacity  (J  kg-1  K-1 ) 

1004 

Air  density  ( kg  irr 3 ) 

1.23 

Oxygen  density  (kgur3) 

1.36 

Nitrogen  density  (kgm-3) 

1.19 

Air  gas  constant  (kg  irr3) 

286.9 

Oxygen  gas  constant  (kg  ur3 ) 

259.8 

Nitrogen  gas  constant  (kg  m-3 ) 

296.8 

Hydrogen  gas  constant  (kg  irr3 ) 

4.1243  x  103 

Oxygen  molar  mass  (kg mol-1 ) 

0.032 

Nitrogen  molar  mass  (kg  mol-1 ) 

0.028 

Hydrogen  molar  mass  (kg  mol-1 ) 

2.016  x  10-5 

Faraday’s  constant 

96485.34 

Universal  gas  constant  (kg  m-3 ) 

8.31451 

pJo(k  +  1)  =  p'„>(k)-“A-(k)  (25) 

9pi,  o 

f)F 

PlUp(k  +  \)  =  plip(k)-a-^(k)  (26) 

dPi,p 

where  a  is  the  learning  factor,  and  p  =  1, . . n. 


3.3.  The  initialization  of  antecedent  parameters 


Since  the  BP  algorithm  tends  to  fall  into  localized  optimization, 
the  initialization  of  the  design  parameter  is  very  important.  In  this 
study,  the  FCM  clustering  method  [29]  is  employed  to  initialize  the 
antecedent  parameter  vectors  U;  and  erf. 


Vi  =  cf 


=  P  x 


EliM- 


X  Xm-Cj 


ElA 


for  i  =  1,  •  • c 


(27) 

(28) 


where  /3(=3)  is  a  given  factor  in  advance,  the  first  set  of  c  data  vectors 
are  obtained  from  the  training  data  set  as  the  initialization  param¬ 
eters.  The  FCM  initialization  method  is  more  likely  to  avoid  local 
optimization  of  the  BP  algorithm. 


with  Hj(k)  = 

Ei=1WiWk)l 

0{k)  =  [/Xi,  ju-2,  Me,  M lXfc,  •  ■  ■,  McXfc,  •  •  •,  Ml^fcn,  •  •  •, 

0=  [P,1.0-"Plc.oP,l.l-"P'c.1-P'l.n"-Pc.n]r 

where  w,- [*(/<)]  is  the  match  degree  of  mth  data  vector 
x(/<)  with  respect  to  A/.  Then  the  consequent  parameters 

[Pi,o'"pc,o  'pii,i"'pc,i--  P,i,n"  pc,n]T  are  identified  by  the  least 
squares  method. 

The  modified  algorithms  significantly  improve  the  accuracy  of 
the  identification  and  ensure  a  more  accurate  model  to  fit  the 
dynamic  characteristics  of  the  SOFC  stack.  It  is  not  necessary  to 
consider  the  SOFC  materials,  structures  and  other  performance 
parameters  in  the  mechanism  model  and  determination  of  these 
factors  is  more  complicated.  Clearly,  the  modified  fuzzy  model  is 
simple  and  convenient  to  build,  as  long  as  the  input-output  test  data 
are  available.  A  multi-variable  model  can  be  readily  established  by 
the  fuzzy  modeling  methods.  In  addition,  after  fuzzy  identification, 
the  fuzzy  model  can  facilitate  the  application  of  an  SOFC  automatic 
control  system. 

4.  Results 

4 A.  Preparation  of  training  data 


3.4.  On-line  adjustment  of  consequent  parameters 

For  a  complex  system  that  operates  over  a  wide  range  of  con¬ 
ditions,  it  is  difficult  to  adopt  an  off-line  model  to  accurately 
approximate  the  dynamic  characteristics  of  the  system.  In  the  pro¬ 
posed  algorithm,  if  the  predictive  precision  of  the  model  cannot 
satisfy  the  need  for  real-time  control,  the  consequent  parameters 
pi  should  be  rectified  on-line. 

Equation  (20)  can  also  be  expressed  as 

E;=iwiMk)] 

C 

=  y>,  W  x  yu(k  + 1 )  =  0(k)  x  9  (29) 

1=1 


In  this  experiment,  the  SOFC  system  is  regarded  as  a  two-input 
and  single-output  system.  The  detailed  parameters  used  for  the 
physical  stack  model  are  given  in  Table  1 .  From  the  physical  SOFC 
model  testing  experiments,  training  data  of  different  air  and  hydro¬ 
gen  flow  speeds  corresponding  to  the  SOFC  stack  temperature  can 
be  obtained,  which  forms  a  data  group  that  is  divided  further  into 
a  training  group  and  a  testing  group.  The  training  group  is  used 
to  identify  the  modified  T-S  fuzzy  model;  and  the  testing  group  is 
used  to  validate  the  modified  T-S  fuzzy  model. 

4.2.  Identification  result  of  modified  T-S  fuzzy  model 

The  training  data  are  classified  into  12  groups  using  the  FCM 
clustering  method,  and  the  results  of  the  FCM  clustering  method  are 
employed  to  initialize  the  antecedent  parameters.  In  the  BP  learning 
algorithm,  the  exponent  m  (m  =  2)  is  a  real  number,  c  (c=  12)  is  the 
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Fig.  6.  Comparison  of  temperature  responses  of  SOFC  stack  generated  by  the  physical 
model  and  the  modified  T-S  model. 


Time(s) 


Fig.  7.  Early  stage  magnification  of  the  temperature  responses  shown  in  Fig.  6. 

number  of  fuzzy  rules,  e  (e  =  0.001 )  is  the  stop  criterion,  the  learning 
factor  is  0.01,  and  the  learning  number  is  100;  as  a  consequence,  an 
off-line  T-S  fuzzy  SOFC  model,  which  consists  of  12  fuzzy  rules, 
is  obtained.  The  comparison  of  stack  temperature  dynamic  curves 
between  the  data-driven  model  (modified  T-S  fuzzy  model)  and  the 
physical  model  (simulation  stack  model)  is  represented  in  Fig.  6, 
with  the  early  stage  shown  in  Fig.  7.  The  temperature  deviation 
curve  is  shown  in  Fig.  8.  Less  than  ±0.6 1<  is  achieved,  suggesting 
that  the  modified  T-S  fuzzy  model  can  approximate  the  dynamic 


Fig.  8.  The  temperature  deviation  between  results  from  the  SOFC  stack  physical 
model  and  the  modified  T-S  model. 


behavior  of  the  physical  SOFC  model  with  high  accuracy  and  provide 
a  satisfactory  control  strategy. 

5.  Conclusions 

In  order  to  facilitate  valid  control-oriented  modeling,  a  nonlin¬ 
ear  modeling  study  of  an  SOFC  stack  using  the  modified  T-S  fuzzy 
model  was  carried  out.  It  has  been  shown  that  the  physical  model, 
which  is  constructed  from  mass-balance  and  energy-balance  equa¬ 
tions  with  MATLAB,  can  be  used  to  describe  the  SOFC  stack  and 
produce  training  data  without  conducting  experiments.  The  tem¬ 
perature  submodel  is  the  center  of  the  physical  model  and  is 
developed  by  enthalpy-balance  equations.  Based  on  the  data-driven 
approach,  the  modified  T-S  fuzzy  model  is  more  attractive  because 
it  avoids  the  use  of  complicated  differential  equations  to  build 
the  SOFC  model  and  quickly  identifies  the  characteristics  of  the 
input-output  relationship.  The  simulation  results  show  that  the 
modified  T-S  model  yields  a  high  accuracy  over  a  relatively  wide 
operating  temperature  range.  If  sufficient  data  from  a  real  system 
are  obtained,  an  optimized  control-oriented  SOFC  stack  model  can 
be  built  automatically  in  the  same  way. 
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